Local observation of antibunching in a trapped Fermi gas 

Torben Miiller/ Bruno Zimmermann,^ Jakob Meineke/ 

Jean-Philippe Brantut,^ Tilman Esslinger,^'[j and Henning Moritz^'^ 

^Institute for Quantum Electronics, ETH Zurich, 8093 Zurich, Switzerland 

^Institut fiir Laser-Physik, Universitat Hamburg, 22761 Hamburg, Germany 

(Dated: 20110413001300) 

Local density fluctuations and density profiles of a Fermi gas are measured in-situ and analyzed. 
In the quantum degenerate regime, the weakly interacting ®Li gas shows a suppression of the den- 
sity fluctuations compared to the non-degenerate case, where atomic shot noise is observed. This 
manifestation of antibunching is a direct result of the Pauli principle and constitutes a local probe 
of quantum degeneracy. We analyze our data using the predictions of the fluctuation-dissipation 
theorem and the local density approximation, demonstrating a fluctuation-based temperature mea- 
surement. 
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PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm, 

A finite-size system in thermodynamic equilibrium 
with its surrounding shows characteristic fluctuations, 
which carry important information about the correlation 
properties of the system. In a classical gas, fluctuations 
of the number of atoms contained in a small sub- volume 
yield a Poisson distribution, reflecting the uncorrelated 
nature of the gas. An intriguing situation arises when 
the thermal de Brogiie wavelength approaches the inter- 
particle separation and the speciflc quantum statistics of 
the constituent particles becomes detectable. For bosons, 
positive density correlations build up, until Bose-Einstein 
condensation occurs, as measured in Hanbury Brown- 
Twiss (HBT) experiments 1.-7- The effect of bunching 
also manifests itself in enhanced density fluctuations in 
real space [§]. In contrast, fermions obey the Pauli princi- 
ple. This gives rise to anti-correlations, which have been 
observed in HBT experiments [HHIS] , and are expected to 
squeeze density fluctuations below the classical shot noise 
limit [T3]. Moreover, for trapped fermions, the reduction 
of fluctuations varies in space, reaching a maximum in 
the dense center of the cloud, which should be accessible 
to a local measurement. 

In this letter we report on high-resolution in-situ mea- 
surements of density fluctuations in an ultracold Fermi 
gas of weakly interacting ^Li atoms. We extract the mean 
and the variance of the density proflle from a number 
of absorption images recorded under the same experi- 
mental conditions. Our measurements show that the 
density fluctuations in the center of the trap are sup- 
pressed for a quantum degenerate gas as compared to 
a non-degenerate gas. We analyze our data using the 
fluctuation-dissipation theorem, which relates the density 
fluctuations of the gas to its isothermal compressibility. 
This allows us to extract the temperature of the system 

[sum US]. 

We first describe the experimental procedure to obtain 
a quantum degenerate gas of about 6 x lO'^ ^Li atoms 
equally populating the two lowest hyperfine states. Fol- 
lowing the method described in [17] , the atoms are loaded 



into an optical dipole trap created by a far off-resonant 
laser with a wavelength of 1064 nm, focused to a 1/e^- 
radius of (22 ± 1) /im ^18;. The cloud is then optically 
moved [12] into a glass cell that provides high optical ac- 
cess, see Fig. fl|a). In the glass cell, forced evaporation is 
performed by reducing the trap power from initially 2 W 
to 4.7 mW. During evaporation a homogeneous magnetic 
field of 300 G is applied to set the s-wave scattering length 
a for inter-state collisions to — SOOag, where ag is the 
Bohr radius. The magnetic field is then ramped to 475 G 
in 150 ms, changing a to —100 oq, and finally the power of 
the trapping beam is increased to 10 mW in 100 ms. Al- 
ternatively, we prepare the lithium gas at temperatures 
above quantum degeneracy. For this, we evaporate to 
50 mW before recompressing to 100 mW, followed by a 
100 ms period of parametric heating. In both cases, the 
cloud is allowed to thermalize for 350 ms before an ab- 
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FIG. 1: (a) Setup for high-resolution imaging of the trapped 
^Li gas. The shadow cast by the atoms held in the dipole trap 
is imaged through the microscope objective and a telephoto 
objective onto an EMCCD chip. The resolution (1/e'^-radius) 
is 1.8 /im [3D] at a wavelength of 671 nm and the magnifi- 
cation is 54. (b) Density distribution (atoms per pixel) of 
the trapped atoms obtained by averaging over 20 realizations. 
The effective pixel size measures 1.2 /xm [20], the maximum 
optical density is ~ 2. The images which we evaluate are 
~ 40 /im wide showing the center of the cloud, which has a 
total length of 500 fim. 



sorption image is taken. The gas is weakly interacting 
since {kpal < 10~^, with kp the Fermi wavevector. 

Our imaging setup is sketched in Fig. [lla). The 
probe hght, resonant with the lowest hyperfine state 
of the |2S'i/2) to I2P3/2) transition, is collected by a 
high-resolution microscope objective and imaged on an 
electron-muhiplying CCD (EMCCD) chip. The atoms 
are illuminated for 8/iS, each atom scattering about 20 
photons on average. Fig. fl|b) shows the average density 
distribution obtained in 20 experiments. 

We now present our procedure for extracting the spa- 
tially resolved variance of the atomic density. The po- 
sition of each pixel in the imaging plane of the camera 
defines a line of sight intersecting with the atomic cloud. 
Correspondingly, each pixel [20] ) having an effective area 
A, determines an observation volume in the atomic cloud 
along this line of sight. At low saturation, the transmis- 
sion t of the probe light through an observation volume 
containing N atoms reads t — e~"'^l-^^ where a is the 
photon absorption cross section. As a consequence, for 
small Gaussian fluctuations of the atom number, the rel- 
ative fluctuations of the transmission coefficient are equal 
to the absolute fluctuations of the optical density and are 
thus directly proportional to the number fluctuations: 



W ^ A^'^'' 



(1) 



where 5t^, (t) and SN"^ are the variance and the mean 
of the transmission coefficient, and the variance of atom 
number, respectively. 

Experimentally, repeated measurements of identically 
prepared clouds provide us with a set of count numbers C 
for each pixel, i.e. each observation volume, correspond- 
ing to a certain number of incoming photons registered 
by the EMCCD. Typically we register ~ 1300 counts, 
corresponding to ~ 130 photons at the position of the 
atoms. We then compute the variance SC^ and mean 
(C). The relative noise of the counts and the relative 
noise of the transmission are related by: 

SC^ _ 2g St^ 



(2) 



Here, g is the gain of the camera for converting photo- 
electrons to counts. The first term is the contribution of 
photon shot noise while the second term is the contribu- 
tion of atomic noise. The factor 2 in the photon shot noise 
term is caused by the electron- multiplying register [21]. 
We extract the contribution of the atoms to the relative 
fluctuations of the counts, by subtracting photon shot 
noise on each pixel according to (pi). This requires the 
value of g (typically ^ 15), which we determine from the 
linear relationship between the variance and the mean of 
the number of counts in a set of repeated measurements. 
The atom number fluctuations are subsequently obtained 
from (fTj). At this stage no division by a reference image 
has been performed, avoiding this source of noise. 



To reduce technical noise adding to these fluctuations, 
we reject images showing the largest deviations of total 
atom number or cloud position [22], which amounts to 
excluding about 30% of the images. The remaining shot- 
to-shot fluctuations of the total atom number SN^^^^ are 
taken into account by further subtracting the quantity 
^Nl^/Nl^{N)^, which is less than <2%oiN ^. Total 
probe intensity variations from shot to shot are below 
0.5%. Applying this algorithm to each pixel of the im- 
ages yields a local measurement of the variance of the 
atom number. 

The mean atom number per pixel is calculated by di- 
viding the mean transmission profile by the mean of ref- 
erence images taken without atoms after each shot, thus 
averaging shot-noise before division. The values for vari- 
ance and mean, obtained by applying the above proce- 
dures, are then averaged along equipotential lines of the 
trap. These lines deviate from horizontal lines (xs-axis) 
in our images by less than half a pixel (0.6 fim). 

Fig. [2] shows the observed variance of the atom num- 
ber plotted against the mean atom number detected on 
a pixel. One set of data was taken for a gas at tem- 
peratures above quantum degeneracy (red squares) and 
another set of data for a quantum degenerate gas (blue 
circles). Above quantum degeneracy, the observed vari- 
ance is found to be proportional to the mean number of 
atoms. The linear behavior confirms that the fluctuations 
originate from atomic shot noise. 
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FIG. 2: Observed variance versus mean of the atom number 
detected on a pixel. Red squares show the data for a non- 
degenerate and blue circles for a quantum degenerate gas. 
The solid red line is a linear fit to the non-degenerate gas, 
yielding a slope of 0.20 ± 0.02. For the data shown, 80 ex- 
periments were performed, 60 for the degenerate case and 20 
for the non-degenerate case. About 30% of the experiments 
were excluded. The error bars shown are estimated from the 
subtraction of photon shot noise which is the dominant con- 
tribution. 

To quantitatively understand the slope of the noise 
curve, which is fitted to be 0.20 ± 0.02, two main ef- 
fects have to be considered. These reduce the observed 
variance and explain the deviation from a slope of one, 
which would be expected for the full shot noise. First, 
the effective size of the pixel is of the order of the res- 



olution of the imaging system. As a consequence, the 
observed noise is the result of a blurring of the signal 
over the neighboring pixels. This effect is also observed 
in [51 [TB], and explains a reduction factor of 0.22 [23]. 
Second, the probe light intensity used for the detection 
is (15 ± 1)% of the saturation intensity. This leads to a 
reduction of the photon absorption cross section due to 
saturation by 0.95 and due to the Doppler-shift by about 
0.9. Together, these effects lead us to expect a slope of 
about 0.19, in good agreement with the observations. 

We now turn to the data taken for the quantum degen- 
erate gas (blue circles in Fig. [2]). At low densities, the 
variance is again found to be proportional to the mean 
density. For increasingly higher densities, we observe a 
departure from the linear behavior and the density fluctu- 
ations are reduced compared to the shot noise limit seen 
for the non-degenerate gas. This is a direct consequence 
of the Pauli principle which determines the properties of 
a quantum degenerate Fermi gas. One can think of the 
Pauli principle as giving rise to an interatomic "repul- 
sion" , which increases the energy cost for large density 
fluctuations. This is similar to the case of bosonic sys- 
tems with strong interparticle interactions, where obser- 
vations have shown a reduction of density fluctuations 
[8] and squeezing of the fluctuations below the shot noise 
limit [Hil]. 
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FIG. 3: Spatially resolved measurement of antibunching. The 
black line shows the mean atom number and the blue circles 
the corresponding variance along the a;i-axis for a degenerate 
gas. The variance is rescaled using the slope fitted in Fig. [2] 
Error bars are estimated from the subtraction of photon shot 
noise, which is the dominant contribution. The dashed line 
shows the variance derived from theory. The shaded region 
indicates the uncertainty originating from uncertainties in the 
trap parameters. The inset shows corresponding data for a 
non-degenerate gas. 

In contrast to previous measurements on antibunching 
[121 113] , we have measured density fluctuations in a spa- 
tially resolved way. For the construction of Fig. [2J we 



have averaged the observed variance for regions of equal 
mean density, whereas Fig. ^ shows the variance (blue 
circles) and the mean (black line) of the atom number as 
a function of the radial position in the trap for a quantum 
degenerate gas. While the variance is proportional to the 
mean in the wings, at low density, we observe a reduc- 
tion of the variance by about 2 dB close to the center, at 
higher density. The inset shows data for a non-degenerate 
gas; in both cases the variance has been rescaled using 
the slope fitted in Fig. [2j The reduction of fluctuations 
is a direct indication of the level of quantum degeneracy 
of the gas. The larger the average occupation of a single 
quantum state, the more the effect of the Pauli princi- 
ple becomes important and fluctuations are consequently 
suppressed. Fig. [3] thus represents a direct measurement 
of the local quantum degeneracy, which is larger in the 
center of the cloud than in the wings. 

To understand this quantitatively, we describe the 
atoms contained in an observation volume in terms of 
the grand-canonical ensemble with a local chemical po- 
tential fixed by assuming local density approximation. 
For a non-interacting gas, the ratio of mean atom num- 
ber and its variance is determined by the fugacity z of 
the system. This leads to the equation 



bN"^ __ \\Jvxii{-z{x\,X2))dxi 

(N) JU-i/2{-z{xi,X2))dX2' 



(3) 



where Li^ is the i-th polylogarithmic function, xi and 
X2 are radial coordinates of the cloud and line-of-sight 
integration is performed along the a;2-axis. 

The dashed line in Fig. [SJis computed using ([s]). For 
the computation we make use of the Gaussian shape of 
the trap, the central fugacity of 13_4 , obtained in an in- 
dependent time-of-flight experiment (see below), and the 
experimental density profile. Our description in terms 
of the grand-canonical potential assuming an ideal gas 
reproduces the experimental data within the error bars 

m- 

We now focus on the interpretation of our results us- 
ing the fluctuation-dissipation theorem. At thermal equi- 
librium, density fluctuations are universally linked to 
the thermodynamic properties of the gas through the 
fluctuation-dissipation theorem, given by 



knT 



d{N) 



5N' 



(4) 



Here T is the temperature of the gas, /i the chemical po- 
tential and fcs the Boltzmann constant. Since the local 
density approximation allows one to assign a local chem- 
ical potential to any position in the trap, it is possible 
to determine the compressibility g ' directly from the 
mean density profiles [TH]. From Q, the ratio of this 
quantity to the measured variance profile of the cloud 
provides a universal temperature measurement [15] . 
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FIG. 4: Fluctuation-based temperature measurement. Vari- 
ance of atom number detected on an efTective pixel versus di- 
mensionless compressibility. The blue circles and red squares 
show the data for the quantum degenerate and the non- 
degenerate case, respectively. The dashed red line is fitted 
to the red squares, giving the temperature according to Q. 



perature using the fluctuation-dissipation theorem |28) . 
In contrast to bosons, fermions cannot exhibit first-order 
long range coherence due to the Pauli principle. How- 
ever, when a Fermi system enters a quantum correlated 
phase, e.g. a superfluid phase, long range even-order cor- 
relations build up [25]. The spatially resolved measure- 
ment of density fluctuations, probing second-order corre- 
lations, is thus a natural tool to study strongly correlated 
Fermi gases. 
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We apply this procedure to our data by computing the 
compressibility, 1 ' = ij [d^j j where we take the 
Gaussian shape of the optical dipole trap into account. 
To avoid the problems of numerically differentiating ex- 
perimental data, we fit the mean density profile with a 
linear combination of the first six even Hermite func- 
tions and use the fitted curve as a measure of the density 
profile in Q [551. Fig. [4] shows the variance of atom 
number plotted against the dimensionless compressibil- 
ity Uq Q ' , where Uq is the trap depth. We observe the 
linear relation described by mn with a slope of ^^ — 
0.27 ± 0.04 for both data sets, the degenerate and the 
non-degenerate. From the physics of evaporative cooling 
it is expected that both slopes are the same [37]. Using 
trap depths derived from the measured laser powers, we 
obtain temperatures of (145±31) nK and (1.10±0.06) fiK 
for the quantum degenerate and the non-degenerate gas, 
respectively. 

To assess the quality of this measurement, we have 
also performed time-of-flight measurements with clouds 
prepared under the same conditions. In this method, 
we determine the temperatures by fitting the measured 
density profiles after free expansion of 1.5 ms (1ms for 
the non-degenerate gas) to the calculated shape of a 
non-interacting gas released from a Gaussian trap. This 
procedure gives us slightly higher temperatures for the 
degenerate and the non-degenerate clouds, which are 
(205 ± 30) nK {T/Tp = 0.34 ± 0.1) and (1.6 ± 0.2) ^K 
(T/Tp = 1.9 ± 0.1), respectively. We attribute the dis- 
crepancy between the two temperature measurements to 
deviations of the trapping beam from a Gaussian shape 
and residual experimental fluctuations from shot to shot, 
which both affect the two methods. 

In conclusion, we have measured density fluctuations 
in a trapped Fermi gas locally, observing antibunching 
in a degenerate gas. We have also determined the tem- 
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